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The microscopic mechanism of circular currents induced in the vicinity of a non-magnetic 

impurity is analyzed for the (time-reversal symmetry breaking) chiral superconducting state 
d(k) = z(k x ± ik y ) which is very likely realized in Sr2Ru04. From the analytic, not self- 
consistent solution of the Bogolyubov-de Gennes equations we find two types of quasiparticle 
states, a bound state and the continuum state. Through impurity scattering the condensate 
transfers angular momentum to the quasiparticle states which generate a circular current. At 
low temperature the continuum part of the quasiparticle spectrum gives the main contribution. 
The non-selfconsistent solution yields a state of finite angular momentum. The comparison with 
the corresponding Ginzburg-Landau theory reveals the existence of a compensating counterflow 
created by the superconducting condensate. The currents and magnetic fields appearing around 
the impurity have possibly been observed by muon spin relaxation measurements. 
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§1. Introduction 



Sr2Ru04 is the first example of copper-free layered perovskite superconductor £F This material 

O '. 

is structurally identical with the layered perovskite La2Cu04, one of the parent compounds of high 

temperature superconductors. In spite of their structural similarity, there are clear differences in the 

electronic properties between the two systems. The cuprates are antiferromagnetic Mott insulators 

in their stoichiometric composition and turn into itinerant electron systems with superconductivity 

only upon doping carriers. On the other hand, the normal state of Sr2Ru04 displays Fermi liquid 

behavior with essentially two-dimensional character and renormalized due to strong correlation 

effects. Because superconductivity appears on the background of strong electron correlation, it 

is unlikely that the pairing channel is conventional s-wave type. It was early suggested that this 
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system might choose spin-triplet p-wave pairing symmetry instead for various reasons.© First 
it resembles the feature of a strongly correlated Fermi liquid like 3 He which is a well-known p- 
wave superfluid. Second, Sr2Ru04 belongs to the Ruddelsen-Popper series Sr n+ iRu n 03 n+ i which 
consists almost exclusively of ferromagnetic compounds, so that the end member Sr2Ru04 may be 
subject to strong ferromagnetic spin fluctuations. 

Indeed there is now overwhelming evidence for the realization of unconventional superconductiv- 
ity in S^RuC^BtJ and the pairing symmetry has most likely the form d(k) = z(k x ±ik y ) analog to 
the (spin-triplet) A-phase of 3 He. The identification of the symmetry is basically possible due to 
two experiments. One is the measurement of the spin susceptibility in the superconducting state 
using 17 0-NMR Knight shift. According to this experiment no change of the spin susceptibility 
appears with the onset of superconductivity indicating equal spin pairing (spin triplet) with the 
moments in the basal plane. EP Note that a strong reduction of the spin susceptibility is common 
to spin singlet superconductors. Furthermore, /iSR zero field relaxation rate experiments show a 
pronounced increase of intrinsic magnetic field in the superconducting phase, suggesting that the 
pairing state breaks time reversal symmetry Tu Both experiments are fully consistent with one 
single state, d(k) = z(k x ± iky), among the possible pairing states. 

The observation of the enhancement of the internal magnetic field is directly related with T- 
violation. In the case of the state d(k) = z(k x ± iky) the Cooper pair has an orbital angular 
momentum parallel to the z-axis of the crystal. This pairing state belongs to the class of "ferro- 
magnetic" T-violating superconducting states according to the symmetry classification by Volovik 
and Gorkov,cP It is natural to expect that the presence of this angular momentum appears in 
magnetic properties. In fact, however, the effect of the angular momentum is invisible in the 
homogeneous superconducting phase due to Meissner screening. It only occurs where the super- 
conducting state is disturbed in some way that screening effects are insufficient. This happens for 
example at the surface of samples and also at domain wall between the phases z(k x + iky) and 
z(k x — iky). Furthermore defects of the crystal lattice, in particular impurities, can also lead to the 
appearance of the unusual magnetic properties of the T-violating state. It is the aim of this paper 
to analyze the effect of scattering by a single impurity on the quasiparticle states and the resulting 
magnetic properties. It is known that spontaneous supercurrents are generated in the vicinity of 
an impurity in such a state .Bt3'lliP This has probably been first investigated by Rainer and Vurio 
in the case of the A-phase of supernuid 3 He.HHP 

Our analysis is based on the solution of the Bogolyubov-de Gennes equations for the single 
impurity problem. In this formulation we can show how quasiparticle states of different angular 
momentum around the impurity are coupled together and give rise to a circular current and mag- 
netic field. This feature can already been observed by solving the problem on a phenomenological 
level with Ginzburg-Landau (GL) theory. Thus before giving a detailed microscopic analysis we 



will introduce the problem by GL theory. 

§2. Impurity problem in the Ginzburg-Landau theory 

We begin our discuss of the superconducting state around an impurity by the discussion within 
the framework of GL theory. This allows us to recognize some of the basic feature of the problem. 
The order parameter of the superconducting state has two complex components so that we can 
write in general the gap function d as 

d(k)=r ?+ d + (k)+r / _d_(k) (2.1) 

with d-t(k) = z(k x ±ik y ). The free energy expansion in these two components has the formli3HS 

/ = ~(\v + \ 2 + l^| 2 ) + \{\^ + k-l 4 ) + 2|r? + |V| 2 
+ |Dr? + | 2 + |Dr?_| 2 + ^{(£»_r/+ )*(£> + 7?_) + (D. V+ )(D + n.)*} + B 2 (2.2) 

where D^ = —iV^/k — A^, D± = D x ± iD y and B = V x A. The free energy density / is 
given in units of B 2 c /Att (B c : thermodynamical critical field), and lengths are measured in units 
of the penetration depth A and B in the unit y2B c = <f>Q / (2tt A£) , and k = A/£. This form of the 
free energy density corresponds to the weak-coupling approach with a single cylindrical symmetric 
Fermi surface. 

The uniform state is immediately obtained by minimizing the free energy with respect to the 
order parameter, 

(77+,^) = 770(1,0) or 770(0,1) (2.3) 

with 1 770 1 2 = 1. Either the state d + (k) or d_(k) are stabilized in the homogeneous phase. The 
gradient terms, however, show us that this needs not to be true in general if the superconducting 
phase is inhomogeneous as we will see immediately. We will in the following assume that the 
homogeneous bulk phase is (t/ + , t/_) = 770(1, 0). 

Let us now add a term to the free energy introducing the effect of the single impurity located at 
r = 0, 



/imp = ^(|77+| 2 + |r ? -| 2 )5(r) (2.4) 



where we use coefficient S to describe the strength of the scattering potential, noting that conven- 
tional impurity scattering is pair breaking for unconventional superconductors, in general, leading 
to a local suppression of the order parameter. We will analyze the behavior of the order parameter 
in the vicinity of the impurity by assuming weak distortion of the order parameter for distances 
sufficiently far from the impurity. Thus we introduce two complex functions «+ (r) and U- (r) with 



r l+ = Vo + u +( r ) an d V- = u -{ r )- Neglecting in addition the effect of a finite vector potential on 
the order parameter we obtain the following coupled GL equation in first order in the impurity 
coupling strength. 

V 2 V 2 

u + (r) + r)lu + *(r) nU+(r) - — 4«_(r) + f] SS(r) = 

(2.5) 
V 2 V 2 . 

where we use V± = V x ± iSJ y . This linear equation system is easily solved by using Fourier 
transformation. If we set tjq = e* x , we obtain the following solution, 

q 2 

roc <_>(! H jj 

e~ lx u + (r) = / dqq ^ ^ 2tv J (qr) = 51 (r) 

•Jo 2 + — o 2 + — -a 4 

(2.6) 
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The important point of this solution is the fact that ?7_-component is induced around the impurity 
and has a phase winding of — 2 x 2ir. Any spatial variation of 77+ drives ??_ which carries a phase 
winding. 

This feature leads now to a spontaneous circular current which we easily find by variation of the 
free energy with respect to the vector potential and the Maxwell equation V x B = j. Since no 
current is flowing in the homogeneous case we express the current j in first order in S, 

j = irj^[a 'Vu + + -(a z + i<j x )Vu-\ + c.c. (2.7) 

where <tq is the 2x2 unit matrix and a^ the Pauli matrices. Also here we neglect the vector 
potential. With the above solution we find the circular current as 
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where a = 2 j=, (3 = 2 -\ — j= and K\(x) is the modified Bessel function of first order. This 

expression shows that the flow direction of the current changes as we move away from the impurity. 



The current generated near the impurity is compensated again farther out. This kind of counter 
flow has been reported also by Rainer and Vurio in the A-phase of 3 He.H3* 

At distances far from the impurity the current decays exponentially, je( r ) = 
—SKa 5 ' 4: ^/3^T^KJzexp(—^/aKr) with a length scale comparable to the coherence length. At short 
distances the expression in eq.( |2.8| ) suggests a 1/r-divergence. This behavior is certainly not ap- 
propriate for short distances. Thus, we need to introduce a cutoff into the integral in eq.( |2.8| ) which 
leads to je(r) oc r for r — ► 0. 

The current and counter current flow in a narrow range around the impurity. The resulting total 
angular momentum of this flow pattern or, equivalently, the magnetic moment vanish due to the 
exact canceling of the two circular currents, 

M = - f d 2 rr x j(r) = (2.9) 

Since the canceling of the current occurs on a length scale shorter than London penetration 
depth (in type II superconductors) the Meissner screening effect does not play an important role 
for the counter flow. Among the two length scales \j K\f]3 and 1/k\/oc, the shorter one describes 
the building up of the magnetic moment due to the order parameter distortion around the impurity 
and the longer one the decay of the r/_-component which due to the phase winding introduces the 
compensating counter flow. Note that this winding does not introduce a topological charge like a 
finite flux of a vortex, since the winding order parameter component does not exist in the bulk. As 
a consequence there can not be a finite magnetic flux associated with the impurity state although 
there is a local magnetic field. 

§3. Bogolyubov-de Gennes Formulation 

We turn now to the description of the impurity problem in the Bogolyubov-de Gennes (BdG) 
formulation in order to understand the origin of the spontaneous currents on the level of quasi- 
particle states. It is possible to perform most of the calculations analytically, if we do not require 
self-consistency of the pair potential, but leave it constant also in the vicinity of the impurity. We 
will discuss the shortcomings in the approximation later. 

The uniform pair potential A(k) = %{a • d(k))a y is given by 

d(k) = A (k x + ik y )z (3.1) 

The BdG equation can be written in the following non-local form 



h (r)u iCT (r) +J2a> J dr' A at(7 <(r,r')v ia/ (r') = EiU if7 (r) 
-h (r)f ifT (r) -J2a> I dr /i* ual {r,v')ui a i(r) = EiU i(T (r) 



(3.2) 



with u\ a {r) and v- 1(J (r) is the particle- and hole-like part of the wavefunction. The part 1iq includes 
besides the kinetic energy also the potential of the non-magnetic impurity located at the origin, 

ho(r)=h kin + U,5(r) (3.3) 

with hkin = — 2m — £f where €f is the Fermi energy. Moreover, U denotes the potential strength 
of the impurity for which we assume the contact type s-wave scattering. 

Because we neglect the spatial dependence of the gap function, in particular the suppression of 
the pair potential around the impurity, we can simplify the equation by Fourier transformation to 
momentum space. 

ikUka + AkVk-a + y Yl Uk ° = eUk ° ( 3 - 4 ) 

k 

-^kVk-a + &*k U ko ~ T7 Yl «fc-<7 = £Vk~a (3-5) 

k 

where A^ = Aocxp(i(/>) (we neglect the gap's dependence on the magnitude of the wave vector |k| 
and take the value of the gap at the Fermi level for simplicity), £fc = k 2 /2m — ep. By introducing 
the variables Io = (U/V) J2k u ka / an d lo = (U/V) J2k v k~a ¥" 0> we can express the wave 
functions as 

(e + £fc)!o- A fc Io' , , 

Uk ° = e 2 - tf - A ^ (3 - 6) 

_ AVq - (e - g fc )I ' 

£ z ~ Kk l ~ A ^ 

Both Iq and I' should be determined self-consistently leading to the equations for the energy, 
tt U ^ (e + £ fc )Io - A fc Io' T U 

Y U ^=yY 

k k 

(3.8) 
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l ° = V l^ V k-° ~ Tl 1^ ,2,2 A2 -^TFZ, 



V k V k ^-^-A 2 - F Z^ e 2_ ?fe 2_ Ao 2 



where for final form we used the fact that the angular integral over the gap function A/, vanishes 
and that close to the Fermi surface the assumption of electron-hole symmetry is satisfied, i.e. the 
normal state density of state is constant. Actually, there is no qualitative change if we include 
particle-hole asymmetry. 

We immediately see that the equations for Iq and I are decoupled and give rise to different 
solutions. The states associated with Io 7^ can be considered as particle-like while the ones for 
I 7^ are hole-like. Within both sectors there are two types of solutions, discrete energy states 
and states belonging to a continuous spectrum (Fig. 2). The former consist of one midgap bound 



state in both sectors and an electron-(hole)-like state above (below) the band top (bottom) which 
correspond to antibound states (for repulsive impurity scattering). In the following we will neglect 
these anti-bound state, since they are very far from the Fermi level and will not have much influence 
on the properties we are interested in. It is useful to separate the discussion of the midgap bound 
states and the states of the continuum. 

3.1 Midgap bound states 

From Fig. 2 we recognize that there is one discrete midgap state for both case Lj 7^ and I 7^ 
with the energies 

^7=f (3.9) 

respectively where c = ttNqU and Nq is the density of states at the Fermi level. Note, that for 
strong scattering (U — ► 00) the bound state energy are both zero and for weak scattering (U — ► 0) 
the bound states are located close to the gap edge ±A(). For the case particle-like bound state 
(Io 7^ 0) we obtain the wave function 

4- -) ( r ) = 7T dkk 2 e ~c 2 k a 2 J °( /gr ) ~ /A ° 2 ° £ ~ 2 h(k F r) 

2tt Jo e- z - £, k - &o VAq 2 - e_ 2 

and for the hole-like bound state (I 7^ 0) 

J o t AU u A °* ;r^J0^ 2I^N A „_,.,„ 



(3.10) 






(3.11) 



where k F is the Fermi wave number, 6 is the angle of position vector r and Jo (At), J±(kr) is the 
Bessel function of 0th and 1st order, respectively. The form of the function fi(kpr) and ^(fc^r) 
for large distances (r >> 1//c_f) is approximatively given by 



v ^0 
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7T / 2 „ 3ti\ 



v^ 1 ^ 



/2(fcFr) "2V^ cos(feFr "T )e "~ " F 



where e = Ag/(\/l + c 2 ) and vp is Fermi velocity. 

The variables Iq and Iq are determined by the normalization condition, 



Ek 2 +0 = i ( 3 - 13 ) 

k 



2 



leading to I 2 = I = A c 3 /(7rN (l + c 2 
The solution in eq.( |3.10| ) and (|3.11| ) are very similar to the bound state discussed in connection 
with magnetic impurities in conventional sup er conductor gl3 03 and was noticed by Buchholtz and 
Zwicknagel.E-T However, the important difference lies in the angular momentum structure of the 
above bound states. We can easily see that the for the electron-like state the particle wave functions 
couples with the hole wave function of an angular momentum reduced by 1. For the hole- like states 
it is just the opposite way around. This property is responsible for the fact that these states can 
carry a circular current. In terms of the wavefunctions u and v the current is expressed as, 

JB(r) = ^- ^[/ (e;K ia (r)Vu ia (r) + (1 - f{ei))v la {r)Vv* i(T (r) - c.c] 

8eAo7rNoC J 2 (k F r) 2 

= ee 7 ==/(e) (3.14) 

my 1 + c 2 r 

where eg is the unit vector of the 6 component, e = Ap/vl + c 2 , and /(e) is the Fermi distribution 
function. We immediately realize that this current disappears for low-temperatures exponentially, 
ex exp(— e/ksT), if U is positive (repulsive). Basically no spontaneous current is caused by the 
bound states if consider the electron states, since the electron state corresponds to an anti-bound 
state due to the repulsive potential. On the other hand, considering the situation from the hole 
point of view the potential is attractive and the hole in a real bound state such that it carries 
circular current. We will discuss this point using the e-h transformation at the end of this section. 

3.2 Contribution of the continuous spectrum 

We turn now to the discussion of the continuous part of the quasiparticle spectrum. Analyzing 
Fig. 2 we immediately see that these states are relevant for the formation of circular currents as 
they are occupied at low temperature. In order to discuss their contribution to the current we have 
to go beyond the discussion we used for the midgap states. A possible formulation can be based 
on the solution of the Lippmann-Schwinger equation for the scattering problem. Let us start with 
the construction of the out-going wave as 

^ + \{r) = <j> k {r) + fdyG k (v,r')\J(r')T 3 ^ + \(r) 



.(r) + G° k (r,0) J^ </> k (0) (3.15) 

K ' ; i-GV0,0)Ur 3 



where t, denotes the i-th Pauli matrix and the function (fr(x) is the solution of the B-dG equation 
without impurity given by 
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™ 1 „ikr 



(3.16) 



for the case Iq 7^ 0, I'o = and 



</>fc0) 



^ (0) fc (r) \ = 1 / A fc 

«(°) fc (r) J ~ y/WE k (E k -Z k ) { E k -£ k 



Jkr 



(3.17) 



for Io = 0, I'o 7^ 0. Note that each momentum k corresponds to two values of the energy, E k 
±\/£ tk 2 + Ao 2 . The Green function entering eq.( [3.15 ) is constructed from these solutions 
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(3.18) 



where r = |r — r |. The Green functions G° k fi(r), G° kt 3(r) and G° kt g(r) can be expressed analyti- 
cally, 
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A;+ 2 k 2 — k- 2 
= 9i,k{r) + ig 2 ,k{r) (3.19) 

where Ji(x) is the first kind Bessel function and Ho' 1 ) (a;) and Ho' 2 '(x) are Hankel functions. More- 
over, the momenta k± are given by 



k± = kF 



\ 



1± 



El ~ A 2 



e F 



k F ± 



El ~ A 2 



v F 



(3.20) 



where vp = kp/m is the Fermi velocity. The two momenta k + and k_ correspond to the more 
particle-like and more hole-like sections of the quasiparticle spectrum. The appearance of these two 
states together is particular feature of the particle-hole mixing in the superconducting state and is 
related with the Andreev reflection. The last of these three Green functions will be important for 
the discussion of the circular current. 
In (3.19 ) we decompose G k> e(r) as gi >k (r) and g2,k{f) as it will be convenient for the later discus- 



sion. It is easy to see that g2,k( r ) is much smaller than g\^ k (r) for the states with \\/E% — Ao| <C ep. 
Since these states will be dominant in the contribution to the currents, in particular, due to the 
enhanced density of states for energies just above the gap, we may safely neglect g2,k(f)- 
The above Green function can now be used to define the T-matrix, 

\Jt 3 

Tie) = s = T t + T 3 t 3 (3.21) 

W l-G°(0,0)Ur 3 { } 

where the two components are given as, 

_ mV 2 p(E k ) 

To - -TTWmW (3 - 22) 

Ts = i + (ME k )uy (3 - 23) 

where p(E k ) = — ImG ^ q(0)- Note that To is purely imaginary, while T3 is real and vanishes in 
the unitary limit. This just reflects the fact that the scattering phase shift is ^ in this limit. 

Let us now rewrite the outgoing wave solution of the Lippmann-Schwinger equation. We restrict 
to the case where Lj = and I'o 7^ since the other case gives finally the same contribution to the 
current. 



V. fc (+)(r 




y/2VEk{Ek — £fc 




(S (r) + S 3 (r))A k 



^2VE k (E k -i k ) \ (5 (r) - S 3 (r))(E k - £ fc ) 



(3.24) 



G° k ,e(r) / (T -T 3 )(E k -£ k )e 



\e 



y/2VE k (E k -^k) \ (T + T 3 )A k e- ie 

with S (r) = G° kfi (r)T + G° K3 (r)T 3 and S 3 (r) = G\ fl {r)T 3 + G\, 3 (r)T . We use now this 
solution to calculate the current according to eq.( |3.14 ). Also in the present case only the angular 



component is finite which we decompose into two part je(r) = j\(r) + J2(r)- The product of the 
third term in eq.( 3.24| ) and its derivative leads to 
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ii(r) =-V E ' '^ (|To| 2 + |r 3 | 2 ) 

fc,.E fc >0 

(3.25) 

^ £ /■ dfcfc 7r 2 iV 2 A 2 (J^fc+Q + J^fc-r)) 2 ^ 2 2 

^L 27r4(£ 2 -A 2 ) r [ll ° l +|J31 j 

where we used that k± « fc_p and g2 <C 51 in eq.( p.l9| ). Next we consider the contribution coming 
from the product of the first and the third term in eq. (|3.24| ), again neglecting the contribution of 

52, 



e r f dkk gi fc(r) m A 
?ni J 2ir hu 



e f dkkg lk (r) J\(kr) 



_e_ y dfcfc ttNqAITq (Ji(fc + r)+Ji(fc_r))Ji(fcr)) _ (3-26) 

mi] Ek >o 2vr 2 . /g2 - A 2 , ^ r 



J fc ^0 



Note that the expression of the current ( |3.25| , 3~2l]| ) the the Fermi distribution function does not 



appear and the only temperature dependence enters only via Ao- In the unitary limit (U — > 00) 
T3 vanishes and the total circular current which comes from the continue parts reduces to rather 
simple form 



. , e f dfcfc A 2 , 



4{Ji(fc+r) + Ji(fc-r)} 2 - J(Ji(fc+r) + J 1 (fc_r))J 1 (fcr) 
lb 2 



(3.27) 



J k 

The main contribution of this integral comes from the states at the gap edge with the wave vector 
close to kp. In this form it is easy to analyze the long-distance behavior of je(r), i.e. r ^> 1/fcj?. 

Mr) « -5^2 [1 + e" 2 ^1 cos 2 (fc F r - ^vr) (3.28) 

where vp is Fermi velocity and both components of the current show the same r-dependence. We 
observe two distinct contributions showing Friedel-type oscillations, one is exponentially fast de- 
caying and the other follows the power-law 1/r 2 (Fig. 3). The latter is a particular result of the fact 
that the current is carried by the quasiparticles belonging to the continuum spectrum. This term 
does not fit well with the solution of the GL theory where we found an exponential decay behavior 
including also the flow of counter currents. While the power-law behavior may be appropriate in an 
intermediate length regime the GL treatment suggests that it should not apply for long distances .t 2 ]) 
Definitely the long-range extension of the circular current should be inhibited by the Meissner effect 
which yields screening counter currents living on the length scale A, the London penetration depth. 
This effect is not included in our discussion here, as it was omitted in the analysis of the GL theory. 
A further more serious omission is, however, the self-consistency of the pair potential A. In the GL 
theory the reaction of the superconductor to the impurity was the admixture of superconducting 
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phase with opposite chirality, i.e. locally the state d(k) = z(k x — ik y ) is admixed to the bulk state 
d(k) = z{k x + iky). 

3.3 Electron-hole transformation 

As mentioned above the repulsive impurity potential problem considered from the electron-hole 
converted point of view corresponds to a problem with attractive potential. This can be easily seen 
by applying electron-hole transformation which means that we replace 

K,h( r ) = ~ n o,e(r), A\(r,r') = -A e *(r,r') (u h ,i, v h>i ) = (v e ,i,u ei i) (3.29) 

Here the subscripts e and h denote the problem in the electron or hole perspective, respectively. It 
is important to notice that the transformation of the gap is actually an inversion of the Cooper pair 
angular momentum, d(k) = z{k x + \k y ) — > z(k x — \k y ). The solution for the so generated attractive 
impurity problem gives the same current as the original repulsive case. 



jh(r) = — 22[f(ei)v* hii(T Vn h i i<r + (1 - /(ei))v/ li)CT Vv^ (T - c.c] 

i,cr 

= je(r) (3.30) 

The main difference between the repulsive and the attractive one is the contribution of the bound 
state solution to the spontaneous current at the zero temperature. In the hole perspective, using 
the bound state solutions ( 3.10) ) and ( 3.11| ), the contribution of the bound state is given as 



8eA 7rN c . ^f 2 (k F r) 2 , . 

jB (r) = -e e === (1-/(6)) . (3.31) 

mvl + c 2 r 

which is non- vanishing. The continuum part compensates for this change, but has basically the 



same structure as eq.( 3.28| ) with slightly different coefficients. Therefore, the sign of the impurity 



potential does not change the spontaneous circular current, but only the view point. 

§4. Discussion 

We have considered a superconducting state whose Cooper pairs have a finite angular momentum, 
a chiral state. Clearly the total angular momentum of the system should be zero in the ground state 
or in thermodynamic equilibrium in the absence of an external magnetic field. Nevertheless, in the 
vicinity of an impurity the intrinsic angular momentum becomes visible in form of a circular current. 
In the Bogolyubov-de Gennes formulation we have seen that these circular currents originate mainly 
from quasiparticle scattering states belonging to the continuous part of the spectrum, while the 
quasi particle bound state at the impurity plays in this respect a minor role. This calculation 
has been done without taking self-consistence into account. Thus the order parameter is uniform 
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and does not get a feedback from the modified quasi particle spectrum. As a consequence the 
quasiparticle states, a hybridization of electron- and hole-like states of different angular momentum 
yield a finite total angular momentum if we do not include the edge of the system. The edge 
introduces chiral quasiparticle modes and a spontaneous surface current .113 In our solution of the 
B-dG problem the impurity potential leads to a shift of the electron concentration and some amount 
of charge is transfered to the edge modes which have zero-energy states. This in turn leads to a 
change of the edge currents and to a compensating change of angular momentum. In our solution 
this is connected with the circular current which decays like r~ 2 . Thus the total angular momentum 
of a finite system is not changed even in the non-self-consistent solution of the B-dG equations. 

On the other hand, the GL formulation leads to a canceling of the angular momentum basically 
within the short length scale £. Thus the feedback of the superconducting condensate leads to a 
strong screening and the impurity and edge problem are decoupled from each other. Thus the above 
connection has to be considered as artificial. This type of screening was recently also discussed in 
detail by Kusama and OhashiEj) in connection with Bloch's theorem. 

While it is difficult to include this screening effect in an analytical way into the B-dG formula- 
tion, we can see this property very conveniently in the GL theory where around the impurity the 
77_ component appears in addition to the bulk phase corresponding to 77+ . This admixed order 
parameter has a phase winding and introduces vorticity yielding a counter flowing current. Since 
this order parameter decays towards the bulk the vorticity does not constitute a real vortex and 
no magnetic flux exists. The r/_ component disappears over a length comparable to the coherence 
length and with it counter currents which within this length completely compensate the angular 
momentum introduced by the impurity scattering states. Overall seen there is no net magnetic 
moment or magnetic flux. Consequently, the circular currents at the impurity are not easy to 
observe. Magnetometers would need a resolution higher than coherence length. One such probe is 
provided by spin polarized muons in zero magnetic field and has recently successfully lead to the 
observation of intrinsic magnetism in the superconducting state of S^RuO^Ef 

Recently, Goryo and Ishikawa have proposed an experiment which has some similarity to the 
impurity problem treated hereHJ Their analysis is based on an topological argument that the chiral 
p-wave state introduces an effective Chern-Simons term into the Ginzburg-Landau theory, which 
couples the vector potential with the scalar potential. A charge inserted into the superconductor 
should generate a circular current similar to a Hall current. There is clear analogies between the 
two studies in various aspects and it is certain worthwhile to investigate the connections more 
carefully. In particular, the argument about the vanishing angular momentum should apply for the 
Hall current of the Chern-Simon term as well in a finite system, since the external magnetic field 
is absent. 

The physics associated with the chiral quasiparticle states around the impurity as well as at 
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the surface can be discussed in terms of Andreev bound states as they were investigated first by 
Buchholtz and ZwicknageliZP for p-wave states and in a related form for d-wave states by HuEiP 
and many others later. In the case of d-wave superconductivity, however, the topological argument 
discussed in Ref.EP does not apply and the zero-energy states at surfaces and impurities are not 
chiral. Nevertheless, they have other important implications £3 

The topological aspect is connected with the presence of an internal angular momentum of the 
Cooper pair of the type p x ± ip y , L z = ±1. Thus, similar physics is expected in other T-violating 
superconducting states of the same or similar topology (angular momentum). To these belong the 
two d-wave states d x 2_ y 2 ±id xy with L z = ±2 and d zx ±id zy with L z = ±1. On the other hand, the 
T-violating state d x 2_ y 2 ± is often discussed in connection with high-temperature superconductors 
does not belong to this class as it does not have an angular momentum. Consequently, we do not 
find a simple circular current around impurities in this case, in contrast to all other case mentioned 
above. The current pattern in this case is more complicated and will be discussed elsewhere. 

§5. Acknowledgement 

We would like to thank A.Furusaki, R. Heeb and T. Tamaribuchi for helpful discussions. This 
work was supported by Grant-in- Aid for Scientific Reserch from Ministry of Education. This work 
is supported by the Grant-in- Aid for Scientific Research (10740169) and (10640341), and one of the 
authors (Y.O) has been supported by a C.O.E. fellowship from the Ministry of Education, Science, 
Sports and Culture of Japan. 



4 

5 

6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 



Y. Maeno Physica C 282-287 (1997) 206. 

K. Ishida, Y. Kitaoka, K, Asayama, S. Ikeda, S. Nishigaki, Y. Maeno, K. Yoshida, and T. Fujita: Phys. Rev. B 

56 (1997) R505. 

A.P. Mackenzie, R. K. W. Haselwimmer, A. W. Tyler, G.G. Lonzarich, Y. Mori, S. Nishizaki, and Y. Maeno: 

Phys. Rev. lett. 80 (1998) 161 

S. Nishigaki, Y. Maeno, S. Farnwe, S. Ikeda and T. Fujita: J. Phys. Soc. Jpn 67 (1998) 560 

G.M. Luke, Y. Fudamoto, K.M. Kojima, M.L Larkin, J. Merrin, B. Nachumi, Y.J. Uemura, Y. Maeno, Z.Q. 

Mao, Y.Mori, H. Nakamura, and M. Sigrist: Nature 394 (1998) 558 

K. Ishida, H. Mukuda, Y. Kitaoka, K. Asayama, Z.Q. Mao, Y. Mori, and Y. Maeno: Nature 396, 658 (1998). 

T. M. Rice and M. Sigrist: J. Phys. Condens. Matter. 7 (1995) L643. 

G.E. Volovik and L.P. Gorkov: Sov. Phys. JETP. 61 (1985) 843 

L.P. Gorkov: Sov. Sci. Rev.A Phys. 9 (1987) 1. 

C.H. Choi and P. Muzikar: Phys. Rev. B. 39 (1989) 9664. 

V.P. Mineev: JETP Lett. 49 (1989) 719. 

D. Rainer and M. Vuorio: J. Phys. C: Solid State Phys. 10 (1977) 3093 

M. Sigrist and K. Ueda: Rev. Mod. Phys. 63 (1991) 239 for a review. 

D.F. Agterberg:. Phys. Rev. B 58 (1998) 14484 

L. Yu: Acta Phys. Sin. 21 (1965) 75. 

H. Shiba: Prog. Theor. Phys. 40 (1968) 435. 



14 



17) L. Buchholtz and G. Zwicknagel: Phys. Rev. B 23 (1981) 5788 

18) Y. Kusama and Y. Ohashi: J. Phys. Soc. Jpn. 68 (1999) 987 

19) M. Matsumoto and M. Sigrist, J. Phys. Soc. Jpn. 68 (1999) 994 

20) J. Goryo and K. Ishikawa: Phys. Lett. A 246 (1998) 549 and preprint |cond-mat 9812412 

21) C.R. Hu, Phys. Rev. Lett. 72 (1994) 1526. 

22) For a review see S. Kashiwaya and Y. Tanaka, preprint. 



15 



Figure caption 

• Fig.l Circular current around the impurity derived from the GL theory. We take the value 
K = 7, and the length is renormalized by the penetration depth A and the current by 4ttSk 2 . 
The expression of ( J2.8Q (solid line) is not valid in short distance and the current should be 
for r — > 0. We omit this in our graph as the length scales are rather different. 

• Fig. 2 Schematic view of the solutions of the B-dG equation. Circles ( diamond ) correspond 
to the continuous ( bound ) states. The solutions in the first and the third (the second and 
the forth ) quadrants are particle ( hole ) like states. 

• Fig. 3 Current around the impurity from eq.( |3.27| ) . We renormalize the energy by Ao, and the 
length by the coherence length £ = vf/Aq, and the current by e/(2m). 
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